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Experiments on the diffusion-limited corrosion of porous copper clusters in thin gap cells contain- 
ing cupric chloride are reported. By carefully comparing corrosion front velocities and concentra- 
Tn . tion profiles obtained by phase-shift interferometry with theoretical predictions, it is demonstrated 

[j^ ' that this process is well-described by a one-dimensional mean-field model for the generic reaction 

A -I- B (static) — > C (inert) with only diffusing reactant (cupric chloride) and one static reactant 
(copper) reacting to produce an inert product (cuprous chloride). The interpretation of the exper- 
iments is aided by a mathematical analysis of the model equations which allows the reaction-order 
r^ , and the transference number of the diffusing species to be inferred. Physical arguments are given 

CIh' to explain the surprising relevance of the one-dimensional mean-field model in spite of the complex 

(fractal) structure of the copper clusters. 
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I. INTRODUCTION 



Q , DifFusion-limited processes are ubiquitous in physics [pj , chemistry M and biology [pi . Reaction-diffusion processes 

have been the subject of intense and continuous interest since the work of Snioluchowski ^^. A crucial feature 

^>~ji of many such processes controlling pattern formation and reaction efficiency is the "reaction front" , a dynamic but 

r*^ ' localized region where reactions are most actively occurring which separates regions rich in the individual reactants. 

Qh, The simplest theoretical model of a reaction front, introduced more than a decade ago by Galfi and Racz H, is the 

"mean-field" model for two initially separated species A and B reacting to produce an inert species C. Since then, the 

,__! ' case of two diffusing reactants A and B has been thoroughly studied analytically Ip-CT and numerically |0-|2^, and 

^ . some predictions of the mean-field model have been checked in experiments p^-pQ . 
^O ' In contrast, the case of only one diffusing reactant A and one static reactant B (confined on a fixed matrix) has 

^^ , not yet been studied experimentally. We show in this paper that the corrosion of a porous solid (B) immersed in a 
^— ^ chemically active fluid suspension (A) can also be described by such a mean- field model. Some analytical |10|3^ and 



i«j numerical n,p2| studies exist for this case as well, but since it is more microscopically complex (for a real porous 

^j. interface) than the case of two diffusing reactants (in a homogeneous medium) an experimental test of the model is 

^\^ \ needed. 

The mean-field model of a planar reaction front for the chemical reaction 



A (diffusing) + B (static) -> C (inert), (1) 

postulates that the concentrations paJX, T) and pb{X^ T) of species A and B, respectively, evolve according to a pair 



Da^-R{pa,Pb) (2) 



f-^ I of coupled partial differential equations |10| , |31| 



dpA ^ „ d'^pA 
dT "^ dX 



C^ ■ -^^ -R{pa,Pb), (3) 

where Da is the diffusion constant for species A and R{pA, Pb) is the reaction rate density. The most frequently used 
initial conditions assume that the reactants are uniformly distributed and completely separated at first, pa{X,0) — 
p'\H{X) and p°g{X,Q) = p'^H{—X), where H{X) is the Heaviside unit step function. Such initial conditions are 
easier to reproduce in experiments than those involving uniformly mixed reactants. There are several assumptions 
behind Eqs. (0)-(||): («) The product C is generated in small enough quantities that its presence does not significantly 
affect the dynamics; {ii) The concentrations are dilute enough that the diffusivities are constant; {in) The fixed matrix 
containing reactant B (static) is porous enough that reactant A can freely diffuse through it; and (w) The reaction rate 
is a function of only the local concentrations and not any fluctuations or many-body effects (which is the "mean-field 
approximation"). It is common to make the mean-field approximation under the assumption R{pA, Pb) = kp^p^, 



but in the interpretation of our experiments we will not assume anything about the form of R{pA, Pb) o, priori since 
the reaction takes place at a solid-liquid interface. Moreover, this interface is highly ramified, and therefore, the 
underlying microscopic dynamics is expected to be more complex than for simple homogeneous kinetics. 

In this paper we carefully test the validity of these assumptions with experiments on a particular porous-solid 
corrosion system: copper clusters corroded by a cupric chloride (CUCI2) electrolyte. The clusters are obtained by 
a thin-gap cell electrodeposition from a CuCl2 electrolyte at fixed current. This process builds a depletion layer of 
CuCl2 ahead of the copper deposit. When the current is switched off, this CuCl2 depletion layer relaxes toward the 
copper cluster bringing Cu^"'' cations which react with copper according to: 

CuCl2 (aq) -f Cu (soUd, red) -> 2 CuCl (solid, white) (4) 

where the cuprous chloride (CuCl) is produced in the form of small (white) crystallites which drop down to the bottom 
of the cell. 

In section II we describe the experimental set-up and the method used to prepare the porous clusters to be corroded. 
In section III we report the experimental evidence that our corrosion system behaves like a ID diffusion-reaction process 
with one static reactant. In section IV, a mathematical analysis is presented which makes quantitative predictions 
based on the experimental data of section II, within the theoretical framework of the mean-field model, Eqs. (||)-(|^). 
In the last section IV, the experimental results are revisited to refine the comparison with the theoretical model and 
to discuss in some detail of its physical limitations. 

II. EXPERIMENTAL METHODS 
A. Apparatus 

The experiments are performed in a thin-gap electrodeposition cell, which is depicted schematically in Fig. IT]. The 
cell consists of an unsupported, aqueous solution of CuCl2 confined to a narrow region of dimensions W = 5cm x L = 
8cm X 5 — 50/^m between two closely spaced, optically flat glass plates (A/4 over 80mmx50mm). Two parallel, 
ultrapurc copper and silver wires (50/im diameter, Goodfellow 99% purity) are inserted between the two glass plates 
to act both as spacers and as electrodes. During the electrodeposition (prior to corrosion) the wires are polarized so 
that the silver wire acts as the cathode and the copper wire as the anode. The solutions of CuCl2 (ACS reagent) 
are prepared from deionized water, carefully cleaned of any trace of dissolved oxygen by bubbling nitrogen through 
it for one hour. The anodic part of the cell (not shown in Fig. nh is filled by a dilute solution of CuCl2 to postpone 
the precipitation of the salt due to saturation effects by dissolution of the anode. The copper electrodeposits are all 
grown at constant current, and the entire experiment is performed at room temperature (« 20° C). 

Digitized color pictures of the copper clusters are obtained by direct imaging of the cluster through a lens, using a 
three-CDD camera coupled with an 8-bit frame grabber from Data Translation driven by the public domain software 
IMAGE [|3| which successively captures three RGB frames and from them reconstructs the color image. 

A phase-shift Mach Zehnder interferometer is used independently to resolve the concentration field, averaged over 
the depth of the cell. A sketch of the interferometer can be found in ref. |Q . The interference patterns are recorded 
through a CCD camera coupled to the same frame grabber |3^ with a 768x 512 pixel resolution. Phase-shift interferom- 
etry offers several significant advantages over traditional interferometry in that it provides an accurate reconstruction 
of the entire concentration field, using a set of successive interference pictures recorded for shifted values of the phase 
difference between two optical wavefronts, and can also be used as an holographic interferometer [ |35| . 

B. Preparation of Copper Clusters by Electrodeposition 

When current flows from the anode to the cathode, charge transfer occurs at the cathode, leading to the reduction 
of copper cations into copper metal according to [ p4| , p6p^ ] : 

Cu2+ + 2e- ^ Cu . 

The actual mechanism of deposition is much more complex than this two electron transfer process since competitive 
reactions involving the solvent species are likely to occur. Nevertheless, in CuCl2 electrolytes, we have observed 
that the formation of cuprous oxide (CU2O) in competition with copper by reduction of Cu^"*" cations is not favored, 
contrary to what is observed in copper sulphate (CUSO4) solutions p8|-p0[|, which can be partly explained by the 
strong adsorption and complexation properties of chloride anions M . This reduction process on the cathode implies a 



local depletion of the copper cations close to the cathode and, therefore, also their replenishment by a global transport 
process, namely diffusion. Although electromigration also contributes to transport, it does not act independently of 



diffusion in regions where electroneutrality is maintained |41|, which means everywhere in the cell outside the 10— lOOA 



thick double layer ||4^,^ . This often misunderstood fact was given a firm theoretical basis by Newman over 30 years 
ago in his asymptotic analysis of the transport equations for a rotating disk electrode p4] , but only recently has it been 
quantitatively verified in experiments (by our group) for the case of constant boundary flux at a fixed cathode ||34| , |36| . 
In summary, the theoretical and experimental evidences indicate that in the absence of convection the concentration 
PA of a dilute, binary electrolyte evolves according to the classical diffusion equation, 

^ = DaV^pa, (5) 

where Da is the "ambipolar diffusion coefficient" for the electrolyte given by a certain weighted average of the diffusion 
constants of the individual ions [^ . 

When the interfacial concentration of metal cation Cu^^ approaches zero, the interface becomes unstable and 
develops into a forest of fine spikes which compete between each other to invade the cell Q,Q. In some cases, a 
"dense-branched" pattern is selected Pq,E^-BG| . This morphology is characterized by a dense array of branches of 
invariant width advancing at constant velocity v through the cell, whose tips delimit a nearly planar front between 
the copper salt electrolyte and the deposit zone. We have shown recently that this growth regime can be modeled via 
a ID diffusion model through the measurement of the copper salt concentration field ahead of the growing deposit by 
interferometry p^JSn] . The experimental concentration field closely fits the "traveling- wave" solution to Eq. (^) , 

Pa{X) ^p°a{i- cxp-^/^-) (6) 

where Ld = Da/v^ and X is the distance to the front edge of the copper deposit, in the direction normal to the front, 
oriented toward the bulk electrolyte. The diffusion length Ld is proportional to P°aI j-, where p\ is the initial bulk 
concentration in copper cations and j is the current density. This diffusion length tends to zero as j I P°a increases, and 
in that limit the concentration profile looks like a step function. Note that p^(X < 0) ~ and pa{X ^ L^) « p\^ 
i.e. the metallic copper deposit leaves behind it a region entirely depleted in copper cations pushing in front of it a 
diffusion layer of constant width extending into the bulk electrolyte. Due to the conservation of copper during the 
deposition process, a linear relation exists between the velocity v of the growth and the interfacial flux of cations J, 
namely vpB = J, where pB is the mean concentration of (metallic) copper in the region of the deposit. 

Using the relation v = DaJL^^ the ratio of the copper concentration in the bulk electrolyte p\ (where it takes 
the form of cupric ions) to that in the region of the deposit Pr (where it mostly takes the metallic form) is easily 
calculated from the basic properties of the electrolyte ||3q , |48| -|5(]|l 

Pb 

where t+ is the transference number |^,Q of the copper cation in a CuCl2 electrolyte. Practically, t+ is a charac- 
teristic of the electrolyte and therefore q will not be a free parameter in our experiments (neither i^ nor q depend 
on the current density j). The closer i+ is to 1, the greater the concentration of copper inside the cluster. In CuCl2 
electrolytes, i+ is expected to be smaller than 0.5, which implies that the copper composition of the deposited zone 
will not go beyond twice the original concentration of CuCl2 in the electrolyte. Therefore, the copper clusters obtained 
by thin gap electrodeposition in CuCl2 are in fact highly porous. 

The large porosity of the deposited copper clusters is of fundamental importance in our subsequent study of the 
corrosion of the copper deposits once the current has been switched off (and the electrodeposition halted) because, as 
a consequence, the cupric ions are able to diffuse freely through the dendrites with approximately their bulk diffusivity 
and then react with a large exposed surface of metallic copper. The low density of the deposit also suggests that 
the product of the corrosion reaction, cuprous chloride (CuCl) crystal, is produced in small enough quantities that 
its presence should not significantly affect the dynamics of the reaction-diffusion process. Therefore, by interrupting 
the current during electrodeposition we can observe a simple reaction-diffusion system with two initially separated 
reactants, copper chloride (A) and metallic copper (B), only one of which is free to diffuse. Since the initial interface 
between the bulk electrolyte and the ramified electrodeposit is planar and the deposit is disordered, it is likely that 
the dynamics of the corrosion process will be effectively "one-dimensional" (ID), in the sense that there might be 
nearly perfect translational symmetry in the two spatial directions {Y and Z) perpendicular to direction of the front 
propagation. Moreover, since the dynamics occurs in three dimensions (as opposed to two for a surface or one for 
a molecular channel), it is also likely that a mean-field, continuum model will be valid, although this may not seem 
obvious a priori in light of the complex geometry of the electrodeposits, which is known to be fractal |51-m]. 



The rest of the paper is devoted to a careful, experimental validation of these hypotheses, showing that our 
system is indeed accurately described by a one-dimensional, mean-field model for the generic chemical reaction, 
A-l-B (static) ^C. We begin in the next section by describing the scaling behavior of the reaction front and accompa- 
nying depletion layer of CuCl2. In the following section, a mathematical analysis of the one-dimensional, mean- field 
model is presented which incorporates the observed scalings and makes quantitative predictions regarding the reaction 
front speed and the concentration evolution. Finally, these predictions are checked with a more detailed analysis of 
the experimental data in the last section, and arguments are given to explain the relevance of the one-dimensional, 
mean-field model for our experimental system. 

III. PRELIMINARY EXPERIMENTAL RESULTS 
A. Temporal Evolution of the Corrosion Front 

At the moment when the current is switched off, the region of the copper deposit is entirely depleted of cupric ions, 
which are thus initially separated from the metallic copper in the deposit. At later times, cupric ions diffuse amidst 
the copper dendrites and react at the metal surfaces, leaving behind cuprous chloride (CuCl) crystallites. In Fig. (a) 
and (b) are shown images of a copper deposit just prior to corrosion and after 30 minutes of corrosion, respectively. 
Note that in Fig. pi (b) the interfacial region between the red copper (the grey color in this picture) and the white 
CuCl is rather flat and thin. 

Focusing on the temporal evolution of this red/white interface, we have observed that, while at first the white layer 
of CuCl appears at the tips of the copper-deposit branches, it gradually becomes fiatter and flatter. As a result, the 
system approaches translational invariance along the Y direction, normal to the growth direction A, thus justifying 
a one-dimensional model for the system involving the single spatial coordinate X (normal to the reaction front). 

By carefully comparing the concentration field of cupric ions obtained by phase-shift interferometry and the 
red/white, Cu/CuCl interface observed on the deposit, the location and extent of the reaction front, where there 
is a significant overlap of metallic copper and cupric ions, can be identified. Following a transient regime (which we 
describe in the last section), it is observed that the reaction front approaches a constant width w ~ T" with a = 0, 



which is consistent with certain mean-field models [ll(j,p_l 31 . Using the theoretical methods pioneered by Galfi and 



Racz m in the case of two diffusing reactants, this scaling was first predicted by Jiang and Ebner |11| using physical 
arguments supported by computer simulations and later by Koza [ pl| using asymptotic analysis. 

Recently, Bazant and Stone |l^ have considered the case of higher-order reactions mA+nB{static) — > C represented 
by the mean-field reaction rate R{pA, Pb) = ^P^Pb ^^^ proved that the scaling exponent for the front width is given 
(uniquely) by the formula 

m — 1 
a = 7r, 7T, 8) 

which holds for any real number to > 1. (The scaling solution does not exist for to, < 1.) In light of this result, 
the experimental observation a = is consistent with the usual one-dimensional, mean-field theory only in the case 
TO = 1. If higher-order reactions were present tti > 1, the theory would predict that the reaction front width increases 
in time (a > 0) although always more slowly than diffusion [a < i). 

The position of the reaction front Xf{T) during the corrosion of a copper deposit (grown from a 0.5M CuCl2 
solution at j = 40 mA/cm^) is plotted in Fig g. Note that, after initial transients have vanished (T > 500 s), the 
reaction front itself "diffuses" with its position given by the scaling, Xf ^ T"' with a — ^, which is is also consistent 



with predictions of the one-dimensional mean-field model ||7|,|ll 31 1. In fact, this diffusive movement of the reaction 
front after long times is a robust feature of all mean-field models for two initially separated species, regardless of the 
reaction orders m,n > 1 or the number of diffusing reactants (one or two), as long as there is no relative advection 
of the two species (due to fiuid flow or some other external forcing) |10| or impermeable membrane to one of the 
reactants [|^. 

B. Temporal Evolution of the Diffusion Layer 

At T = when the current is interrupted, the reactants Cu and CuCb are completely separated, since the concen- 
tration of CuCl2 is negligibly small in the immediate vicinity of the metallic Cu electrodeposit. During the subsequent 
corrosion process the concentration of CuCl2 remains very small in the reaction front, which leads to the modification 



of the initial depletion layer of CuCl2 (produced by the electrodeposition process) into a region where the concentra- 
tion smoothly interpolates to the value of the bulk solution far behind the front. The term "diffusion layer" is used to 
describe this region because it is characterized by the transport of fresh CuCl2 by diffusion from the bulk, relatively 
unaffected by chemical reactions due to the negligible (or vanishing) concentration of metallic Cu remaining behind 
the reaction front. 

The temporal evolution of the diffusion layer is revealed by precise interferometric measurements of the concentration 
profile of CuCl2. In Fig. ^ are shown three contour plots of the concentration field computed from the integrated 
index along the depth of the cell. Since the experiments are performed in thin-gap cells (50/im) and the depletion 
layer spreads over distances larger than this gap, it is safely assumed that the concentration of CuCl2 does not change 
appreciably along the Z direction (parallel to the laser beam Q). If Fig. ^, the shadow of the Cu/CuCl cluster is 
also clearly seen. A close inspection of the panels (b) and (c), which correspond to eroded clusters, reveals that in the 
zone of the copper cluster where CuCl2 has diffused (recognizable where the leftmost isoconcentration contours have 
moved through the cluster), the cluster has been broken down in smaller crystallites, which, as indicated by their 
color in Fig. ||, are made of CuCl. 

Typical experimental concentration profiles of CuCl2 measured at different times (averaged along the Y direction, 
normal to the growth direction) are shown in Fig. 0. The shape of these concentration profiles is discussed in the 
next two sections, but here we focus on the scaling of the width Wc;(T) of the diffusion layer (defined as the region 
of non- negligible gradients). Fig. g shows that at long times (T > 500 s), the diffusion layer approaches a self-similar 
structure, with the diffusive flux entering the reaction front obeying the scaling law, J 4^ ex {d p aI dXy\x=Xt '^ T~ , 
and that, therefore, the width of the diffusion layer has the familiar scaling |7|,P, |ll| , |31| Wd ^ T^ with & ~ \, which is 
another robust feature of the mean-field models |jlO] . 

A physical argument based on mass conservation between the diffusion layer and reaction front \^^^ can be used 
to predict the scaling of the reaction rate (per unit volume) in the front R ~ T^^ from the preceding experimental 
observations. The total reaction rate in the front (per unit area) scales as wR ^ j^a-p^ ^^^^]^ \i^[^ Q^\yi of cupric 
ions due to reactions must balance the diffusive flux entering the front Jd '^ T"*, which yields the scaling relation, 
f3 = a + S — 0+^ — ^. It is important to point out, however, that while a — and S = ^ are the results 
of direct experimental observations, the scaling exponent /? = ^ is only inferred by a physical argument, based on 
the assumption that chemical reactions are negligible in the diffusion layer. Although this assumption has been 
checked numerically and analytically for various mean-field models, the reaction rate is not directly measured in our 
experiments. 

In the general case R{pA, Pb) — kp^p^j^ mentioned above, it can be shown jlOJ that /3 is given (uniquely) by 

P=^, (9) 

771-1-1 

SO that once again tti = 1 is suggested by the inferred value /3 = i However, given that the experimental system 
has complex fractal structure and three-dimensional transport in the reaction front, it is not obvious a priori that 
R{pA: Pb) = f^PApB is ^ reasonable approximation within a spatially averaged one-dimensional model. Instead, we will 
make no ad hoc assumptions about the functional form of the reaction rate R{pA, Pb) and then explore consequences 
of only our direct experimental observations, a = and S = a = ^, within the framework of a one-dimensional 
mean-field model. 

IV. THEORETICAL PREDICTIONS OF THE MEAN-FIELD MODEL 
A. Dimensionless Model Equations 

The model equations have a dimensionless form involving only the parameter, q = p'^/pg, defined in Eq. (|^), 

f^='qria,b), (11) 



with boundary and initial conditions 



a(-oo,t) = 0, a(oo,t) = l (12) 

5(-oo,t) = l, 6(oo,t) = 0, (13) 

aix, 0) = H{x), b{x, 0) = H{-x) (14) 



where 



a^P±^ b^P^, (15) 






P^ : ^ - ^ y (^^^<, ) • (17) 

These initial conditions are closest to the actual ones used in the experiments when the copper deposit is grown at 
large current, which corresponds to small L^ in Eq. (0). The initial-boundary- value problem of Eqs. (p^)-(p^ involves 
an idealized, infinite system possessing no natural length or time scale, and, therefore, it is expected that asymptotic 
similarity solutions exist in which distance and time appear coupled by power-law scalings p5[ . The experimental 
system, on the other hand, possesses several relevant length scales, but they turn out not to affect the evolution of 
the reaction front, at least for some range of times. For example, the spatial scales of the copper deposit, such as the 
typical dendrite spacing and dendrite width, surely affect the dynamics at early times since these length scales are of 
the same order as the diffusion length Ld [pq] , but it is observed that during corrosion the system quickly approaches 
planar symmetry, averaged across scales much larger than individual dendrites. Likewise the length scale of the gap 
spacing is not expected to greatly influence the corrosion dynamics because vertical (buoyancy-driven) convection, 
which has been observed during the growth phase |p9] is suppressed in 50 /im depth cells [p4 57|. However, the 



settling of the reaction product, CuCl crystallites, could have some effect on the front dynamics at this scale. Finally, 
the largest length scales, namely the distances from the outer edge of the deposit to the two electrodes, also should 
not affect corrosion dynamics until the reaction front gets close to the cathode and/or the diffusion layer approaches 
the anode. Therefore, during intermediate times, after three-dimensional transient effects have subsided but before 
the system size begins to matter, the corrosion dynamics should be well described by a self-similar solution to the 
one-dimensional mean- field equations. 

B. The Diffusion Layer 

Motivated by these arguments and the experimental data, we consider the transformation 

a(x, t) = i(C, t), b{x, t) = B(C, t), where C = "" ~ ^^^\ (18) 

for the concentration of CuCl2 in the diffusion layer (defined by C > 0), and seek an asymptotic similarity solution, 
A{C,,t) ^ A{C,) and B{C,.t) ^ B{C,) with power-law expressions for Xf{t) and Wd{t). The experimental observations 
discussed in the previous section support the scaling law Wd ~ i^ " for the diffusion layer width and a similar diffusive 
scaling law for the reaction front position Xf ^ t^". Therefore, we make the definitions 

Xf{t)^~2vVt , (19) 

C = ^ + .. (20) 

where v{qf' is an effective diffusion constant for the reaction front to be determined during the analysis. 
Substituting these expressions into Eq. (lOh, we have, 

which simply amounts to a change of variables from {x, t) to (C, t). We now look for an asymptotic similarity solution 
by assuming that the time derivative vanishes relative to the diffusive term, 



dA 
lim t ^- = 0, for C > fixed, (22) 

t— ►oo Ot 



which has been called the "quasistationary approximation" in the physics literature P,15 31 58]. This is not really an 



"approximation" but rather is an exact asymptotic property of a certain class of self-similar solutions to Eqs. ( 10)-(|lJ 
which happen to accurately fit the experimental data, as we will show in the next section. At this point it is common 
to assume that the reaction term also vanishes relative to the diffusion term, 

lim t ■ r{A, B) = 0, for C > fixed, (23) 

t — >oo 

and that the concentration of the non-diffusing species also vanishes in the diffusion layer, i.e. where the reaction 
front has already passed, 

lim S(C;t) = 0, forC>Ofixed. (24) 

t — >oo 

Note that -B(C, t) — Qior C,> v a± all times due to the initial condition 6(a;, 0) = H{—x) and the fact that this reactant 
does not diffuse. The limits in Eqs. (|2^) and ( p4| ) have previously been taken as ad hoc assumptions [gl|, but it can 
be shown that they are in fact necessary consequences of the quasistationarity [Eol 



Using Eqs. (p2|)-(p3[) and passing to the limit t — > cx) with C, > Q fixed in Eq. (21) yields an ordinary differential 
equation for the asymptotic diffusion-layer concentration ^(C), 

A" + 2{Q-v)A' = Q. (25) 

The solution to this equation subject to the boundary condition A{oo) = 1 can be written in terms of error func- 
tions p^ , 

. , >N , , . N erf(C — v) + eviiv) , , 

^(C) = A + (i-^o)^\:^^^, (26) 

where Ao = ^(0) is a constant to be determined by asymptotic matching with the reaction front as C ~^ 0. The 
function A(C) is shown in Fig. R for different values of v. The slope of A(C) at C = given by 

V7r(l -l-erf(i^)) 

is the (dimensionless) diffusive flux into the reaction front. 

On the length scale Wd{t) ^ t^^"^ appropriate for the diffusion layer, the self-similar asymptotic concentration fields 
just derived appear not to be differentiable at (^ = 0, 

a{x, t) - A{C)H{0, b{x, t) - ff (-C), as i ^ oo with C 7^ fixed, (28) 

but, as we have already observed experimentally, that is only because in reaction front (at Q = Q) the concentrations 
are smoothly interpolated across these apparent discontinuities on a much smaller length scale w ^ t" = o{Wd) since 
a < 5. In mathematical terms, the asymptotic approximations in Eq. (Bq) are not uniformly valid for all (x, t) as 
t — > 00, but rather are valid only for C, ^Q, i.e. \fl = 0{\x + 2i'\/t\). 

C. The Reaction Front 

We now explore the consequences of the experimental results a = and S ^ a = ^ within the present mathematical 
model. Although the physical arguments made above for the lack of a natural length scale are much more tenuous in 
the reaction front because the observed front width (about 0.2 mm) is comparable to the average dendrite thickness 
(0.1 mm) and spacing (0.4 mm) as well as the gap (0.05 mm), the nearly perfect planar symmetry of corrosion process 
leads us to nevertheless seek another asymptotic similarity solution to the one-dimensional, mean-field equations in 
the vicinity of the reaction front, x — Xf{t) = 0{1). The predictions of the model will be carefully tested against the 
experimental data in the next section. 

Since a = and a ^ ^, we consider the transformation 

a{x,t)=t-''A{T],t), b{x,t)=B{r],t), (29) 



where 77 is a new similarity variable for the reaction front defined by 

r] = x + 2vt^l'^ = 2t^/^C. (30) 

The exponent 7 > is introduced to allow for the possibility that a{x,t) — > in the reaction front, which is suggested 
by the result r{a, b) ^ t^^ with /3 = ^ inferred earlier from the experimental data. In contrast, no such prefactor 
multiplies B{r], t) in the reaction front since b{x, t) must remain finite there in order to interpolate between the limiting 
values of and 1, respectively, behind and ahead of the front. 
Making these transformations in Eq. ( pO| ) yields 

dA u?dA , r d^A 



"^*"'-^= a;;^^''"('"'-^'^)- ^^^) 



dt drj drj 

As before, we explore the possibility of self-similar quasistationarity in the reaction front: A{i], t) ~ A{'q) and B{ri^ t) ~ 
B{ri) as t — » CX3 with \rj\ < cx3 fixed. The consequence of the quasistationarity assumption in Eq. (pTl) is 

A"{f])^ lira f r(t-"<A{r/),B{T])) , for fixed r/ (32) 

Since A"{ri) — cannot satisfy the boundary condition ^(—00) = (except in the trivial case A{r]) = 0), the limit 
on the right-hand side of Eq. (P2) must be nonzero (and finite), which is possible only if r{a, h) is linear in a, i.e. 

r(a,6)=a/(6), (33) 

for some function /(6). Therefore, the experimental facts w{t) ~ i" and Xf{t) ^ t^''^ are consistent with the one- 
dimensional, mean-field model only if the reaction rate is first order in the diffusing species. 

Next we make the same transformation in Eq. (0) and replace the reaction term with Eq. (p3) to obtain 



f +--"'! = -.'- W(«). (34) 

By inspection, quasistationarity is possible only if 7 = i, which would imply r{a, b) ^ t~^f'^A{'q)f {B{r])). Therefore, 
we conclude j3 = ^ once again, and the physical argument given in the previous section is found to have sound 
mathematical justification. 

With these results we arrive at a third-order system of nonlinear ordinary differential equations for the concentration 
fields in the reaction front, 

A!' - Af[B) = 0, (35) 

vB' + qAf{B) = 0. (36) 

These equations may be combined to eliminate the reaction term and integrated once using the boundary conditions 
ahead of the front, A{—oo) = and B{—oo) = 1 to obtain 

qA! = u{l - B). (37) 

Before proceeding with another integration, however, a third boundary condition is needed, which comes from asymp- 
totic matching with the diffusion layer. 

D. Asymptotic Matching 

In mathematical terms, our equations possess an "internal boundary layer" ||60(| . The reaction front, defined by 
|a; -I- 2v\/t\ = 0(1), acts as the "inner region", while the diffusion layer, defined by ^fi = 0{\x + 2i>^/t\), acts as the 
"outer region". For consistency, the "inner limit" (C -^0) of the outer approximation, Eq. (p8|), must match the 
"outer limit" (77 — > 00) of the inner approximation, Eq. (|29|). We have shown that 7 > is required to describe the 
experimental data, which means that a{x, t) approaches zero uniformly in the reaction front. Therefore, by matching 
at zeroth order we obtain ^'(0) = Aq = 0, but this does not provide the missing boundary condition for the reaction 
front. At the next (linear) order we have 



dA ac, ^^ A'io 
^~l^l^li-'^' as i -> oo, with H<oo fixed ^^^^ 



da I If |§ -- ^, as t -> oo with < C < oo fixed 



^/t dri dx yt 

and by matchin g w e conclude A'{(X)) — Ai, where -4i(j^) = A'{0)/2 can be expressed in terms of ^{q) using Eq. (|27 
In light of Eq. (M), the matching condition for b(x,t) is S(oo) = 0. 

The matching conditions allow us to derive an exact expression for j/(g) and hence the asymptotic front position 
Xf{t) — 2iyy/i. Taking the limit 77 ^ oo in Eq. ( p7| ) using A'{oo) = Ai and B{qo) = 0, we obtain qAi{v) = v. By 
substituting Ai{v) from Eq. (|27| ) we obtain the desired expression for v{q), 

V = F-^{q), where F{x) ee ^/¥xe''' [1 + erf(a;)] , (39) 

which has also been derived by Koza [pll. The relation q — F(y) is plotted in Fig. pi and will be used in the next 
section to estimate q from the experimentally measured value of v. 

With these results, we are led to a second-order, nonlinear boundary-value problem for the reaction front concen- 
tration of the diffusing species: 

A' =A!{\-AIAx), yl(-oo)=0, y^'(oo)=A. (40) 

Note that Eq. (Eo) is invariant to translation j] ^^ r] + rjo, where rjo is an undetermined constant depending on the 
initial conditions that precisely defines the location of the reaction front {e.g. as the point of maximal reaction rate). 
Since it is difficult to accurately measure the reaction- front concentration fields in our experiments, we stop here 
and refer the reader to the article of Bazant and Stone |lO| for the integration of this boundary-value problem and 
other analytic results in the case /(6) = 6™, to > 1. 

V. EXPERIMENTAL TEST OF THE THEORETICAL MODEL 

A. Check of the exact asymptotic predictions 

In section III we showed that as corrosion proceeds the reaction front moves with the time as Xf{T) ~ T'^/'^ 
and does not spread {w{T) ~ T" with a = 0) and the width Wd of the depletion layer increases with the time 
as Wd(T) ~ T^/^. In section IV we showed that these observations are consistent with the predictions of a one- 
dimensional A-l-B (static) -^ C (inert) mean- field model with a reaction rate that is first order in the diffusing species 
A. By solving the mean-field equations, we derived not only the scaling exponents for Xf{T) and W{T) but also 
the prefactors and the exact asymptotic shape of the concentration profile of the diffusing reactant as a function of 

the reduced coordinate C — ,I^L . In this section, we quantitatively test these theoretical predictions against the 

experimental results. 



In dimensional units Eq. (|l^) reads: 



1. Movement of the front 



-Xf ^ 2iy{q)VDf. (41) 



Therefore, from a log-log plot of Xf as a function of T one gets the value of i^, and q can then be deduced from 
Eq. (|39|). In our experimental system, q is linearly related to a characteristic property of the electrolyte, namely the 
transference number of the cation, through q ~ 1 ~ t+. To derive the values of q and t^ from Eq. (pl|), we need 
an accurate value of the diffusion coefficient of the electrolyte. D is likely to depend on the concentration of CuCl2 , 
but to our knowledge, has not been tabulated for CuCl2. Hereafter, we use the value D — (1.0 ± 0.1)10^^cm^-s~^, 
determined independently by our interferometric technique. 

The two sets of experimental data in of Fig. y give ^v^pD ~ (1.7 ±0.1) 10~'^cm-s~^/^, therefore v = 0.27 ±0.02 and 
t+ = 0.33 ± 0.05 from Eq. (|39|), q = F(y). Note that t+ « 0.3 (for a 0.5 mol.l"-^ electrolyte) is quite consistent with 
the corresponding value at infinite dilution t°^ = 0.4 since i+ is likely to be a decreasing function of the concentration 
p9[ . Although we have not directly measured the transference number i+ of the Cu^+ cation, its reasonable value 
just inferred from the observed front speed via Eq. ( |39|) constitutes a successful prediction of the one-dimensional 
mean-field model. 



2. Width of the depletion zone and whole concentration profile 

In this section, we analyze the experiments performed with a higher electrolyte concentration, namely 1.0 mol.l^^ 
CuCl2 . The concentration profile in the laboratory frame can be written in dimensional units using Eq. ( |26| ) and the 
definition of C : 

ejV^/2Vm+^rii^^ (42) 

1 + CTl[l>) 

Note that a{X,T) is used in the experimental parts to denote pa{X,T)/p'\. A characteristic feature of these profiles 
(and the experimental data in Fig. 0) is that they exhibit a fixed point with ordinate: 

aiX = 0,T)^-^^ . (43) 

1 + ert(z^) 

Since a(0, T) depends only on q, a value of q can be deduced from Fig. 0, which shows the concentration profiles 
during the corrosion of a copper deposit obtained by electrodeposition from a 1.0 mol.l^^ CuCl2 solution. We find 
a{X = 0,T) = 0.25 ± 0.01 which imphes v = 0.30 ± 0.01. From Eq. (|3|) the mean-field model would predict 
q — 0.79 ± 0.06. As expected, the inferred value of the transference number, t^ = 1 — q = 0.21 ± 0.06, for this 1.0 
mol.l"^ CuCl2 solution is lower than the value of 0.33 ± 0.05 at 0.5 mol.l"^ computed above. This value is somewhat 
smaller than expected based on concentration effects (see below). Note that we have not directly measured the ratio 
q — p°a/ p% or the transference number t^ in the experiments described in this paper, but the value oi q = 0.79 
just obtained from Eq. ( [43|) is necessary for comparison with the mean-field model (without any other adjustable 
parameters). Therefore, we will use q — 0.79 in the following analysis of the experimental runs in 1 mol.l^^ CuCl2 
electrolyte. 

^From Eq. ( p6| ) the width Wd of the diffusion layer (with dimensions) is given by: 

From an experimental point of view, it is simpler to measure a{X,T) at A" = rather than at A" = Xf{T), so we 
consider the temporal evolution of the gradient of a{X,T) at A' = 0. From Eq. (p2|) we obtain: 



da{X,T) 



dX 



1 



x=o V^DT{l + eTi{v)) 



(45) 



and Wd{T) — exp(z^^)/ dxa{X, T)|y „. Figure ^ shows the quantitative agreement between the experimental values 
of dxa\x=o ^'^'^ ^^^ function of Eq. ( p5[ ) plotted for D ~ 10~^cm^-s~^ and q — 0.79. Note that D and q are deduced 
from previous analysis and are not adjustable parameters. 

Continuing our quantitative analysis of the experimental concentration field, we plot in Fig. g the asymptotic shape 



of the concentration profile. To determine a{Q from a{X,T), we compute C, using ^ — — ^= + v{q), with q — 0.79 



jiiit;. _LU Lit^Leiiiiiiit; uyt^) iiuiii ciyyv ^ ^ y, we cuiiipuie ^ Ubiiig t^ — — 

and D = lO^^cm^.s^^ and adjust the origin of the abscissa to the initial front of copper position, to ensure that 
A{( = 0,T) — for all T. For comparison we also show in the same plot the theoretically predicted function A{() 



function computed from Eq. (|26| ) and ( |39| ) with q = 0.79. 

To focus on the region of the reaction front, the experimental data is plotted in fig nO according to the linearized 
version of Eq. ( |2q ) 

X -Xf 2^""^ X-Xf 



aVDT = A\0)^—^= \ (46) 

2 v'"'!! + ert(i^)) 2 

Since {X — Xf)/2 is proportional to the reaction-front similarity variable rj in Eq. (pO|), the mean-field model would 
predict a collapse of this data to a single curve given by the solution of Eq. (Eol) . 

Unfortunately the noise in the experimental data washes out the exact concentration profiles in the reaction front 
on this scale, but it is clear that the width of the reaction front has the asymptotic scaling w ^ t" with a = 0. 
Moreover, the asymptotic shape of the concentration distribution is quite consistent with the solutions to Eq. ( |40| ) 
given in Ref. |10|. Note that the decay of slope the reaction-front concentration A'irj) toward its limiting in the 
diffusion layer Ai in Fig. nS appears to be quite fast. If this decay were exponential rather than a (much slower) 
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power law, then according to the mean-field model ||T^ the reaction rate would have to be first order in the static 
reactant m = 1, i.e. fjb) = b or r{a, b) = ab, but it is impossible to reach this conclusion definitively from our data. 

As shown in Figs || and |^, all of the measured concentration profiles collapse to the single asymptotic curve 
predicted for q — 0.79 over the whole length scales investigated in the experiment. This quantitative agreement 
between the experimental and theoretical concentration profiles of the diffusing reactant independent of the length 
scale strongly support our modeling of this corrosion experiment with a one-dimensional A + B (static) — > C (inert) 
mean-field model. 

B. The transient 

The A -I- B — > C mean-field model with two diffu sing reactants exhibits many surprising and nontrivial behaviors 
at short times (see [Q and references therein, [Q, fGlJ). In this case, some microscopic parameters like the reaction 
constant(s) can be determined from these short time behaviors. In particular, at a time inversely proportional to the 
microscopic reaction constant, the global reaction rate switches from an initial i^/^ increase to a subsequent t^^/^ 
decrease IQ- Moreover, in the reversible A -I- B ^ C system, a crossover between irreversible and reversible regimes 
can be observed at long times pW and the value of the backward reaction constant can be inferred from the crossover 



time 1 26 . 



In the present case of one static reactant, it is also possible to express the transient decay to the asymptotic solution 
in terms of the reaction orders m and n for the one-dimensional mean-field model pOf . In our experiment, however, 
the transient behavior is determined by a superposition of different mechanisms since our system is not really one- 
dimensional or homogeneous. We now show that the transient behavior appears to be governed by two-dimensional 
geometric effects that hide the kinetic features by analyzing the detail the experimental runs. 

Looking at Fig. Q(a), note that the concentration field is not one-dimensional at the early stages of the corrosion 
experiment: the isoconcentration lines closely follow the jagged outline of the deposit in the region near the tips. 
The amplitude G of the modulation of the leftmost isoconcentration line (the closest to the copper cluster) is about 
0.4 mm. This system clearly cannot be viewed as one-dimensional until the front has traveled at least a distance on 
the order of G. In Fig. |^, note that the time of the transient regime (before the asymptotic i^/^ behavior sets in) 
closely corresponds to the time needed for the front to move across a distance G ^ OAmm. (This two-dimensional 
geometric effect also may explain why the initial movement of the front is slower than the asymptotic behavior, as 
shown in Fig. ||.) 

To further support this hypothesis, we now study the relaxation dynamics of the concentration field. In Fig. yjj- 
(a), is plotted the isoconcentration line corresponding to a = Pa/Pa = 0-lj just after the current has been switched 
off. This line is not continuous, because the concentration field cannot be extracted by interferometry in the zones 
containing the deposit. This line defines a function X(Y) roughly periodic, of amplitude G{T) and period A ~ 1mm. 
It is reasonable to expect that the characteristic time for the relaxation of this modulation of the concentration field 
toward a flat two dimensional profile is the time r/ needed for the front of copper to move from its starting position 
{Xf(T = 0)) on the length scale G{T = 0) = 0.3 mm= 2v^/DFf, which yields the estimate r/ = G{^Yl^v'^D - 250s. 
Moreover, in light of the analysis of Krug |62l described below, it is also reasonable to expect that the functional form 
of the decay will be exponential. 

In Fig ^l]-(b), we plot log(G(r)/G(0)) as a function of the dimensionless time T/rf. The relaxation is well fitted 
by an exponential function, with a characteristic time close to r^, which supports our hypothesis. Therefore in our 
experiments, the transient behavior is directly linked to the relaxation of the initial two dimensional concentration 
field towards a Y-invariant profile and cannot provide information on the kinetics independently. 

C. Physical relevance of the one-dimensional mean- field model 

In the previous sections, we have demonstrated the quantitative agreement between the behavior of our thin-gap 
corrosion system and various predictions of a one-dimensional mean field model. This agreement is not obvious a 
priori, and therefore we close in this section by giving physical arguments to explain this surprising fact. 

1, No inhibition of diffusion or reaction by CuCl 

In Fig. U we see that the product of the reaction does not seem to disturb the concentration field of the diffusing 
reactant A. To understand this fact, we consider the volume occupied by the product CuCl in the cell. We know 
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from Eq. (Q) that the mean concentration of copper before the dissolution is w 2p^. We deduce from Eq. (Q) 
that if CuCl does not diffuse (which is verified in our experiments), the mean concentration of CuCl is twice the 
initial concentration of copper, which is approximately four times the initial concentration of CuCl2 in the bulk, i.e. 
2 mol.l""'^. Since the density and molecular weight of CuCl are 3.38 g.cm~^ and 99 g.mol"'^ respectively, the volume 
occupied by the solid CuCl after the dissolution is roughly 5% of the total local volume. Therefore, the small crystals 
of CuCl do not significantly alter the volume free for the diffusion of CuCl2. Moreover, because the CuCl crystalites 
do not adhere to the copper metal branches and fall to the bottom of the thin gap cell, the surface of the copper 
cluster is constantly renewed and "ready" for corrosion by CuCl2. 

2. Stable front, asymptotically one- dimensional 

The fact that the dissolution process builds a stable (flat) interface can be understood by considering that diffusion- 
limited corrosion is the "time-reversed" process of diffusion-limited aggregation and that the fluctuations of the inter- 
face decay rather than grow to reach a stable flat front asymptotically. Krug |6^] showed that periodic perturbations 
of a flat front of wavelength A in the direction perpendicular to the direction of motion of the interface would decay 
with a characteristic time r = X/v. The stability of the corrosion front can therefore be qualitatively understood with 
the following argument: the electrolyte most easily reaches the most exposed or least screened parts of the copper 
deposit. These bulges are dissolved first, and the interface is smoothed. 

3. Relevance of ID approximation of the concentration field 

In the long-time asymptotic regime, the modulation of the initial concentration of reactant A (CUCI2) relaxes toward 
a flat concentration profile along the direction Y whose shape is given by Eq. (Eq). However, the concentration of 
the static reactant B (Cu), as well as the concentration of the product C (CuCl) keep a periodic shape along the Y 
direction, which somehow does not alter the one-dimensional asymptotic solution. The largest characteristic length of 
the deposit in the direction parallel to the front (Y) is the mean distance A between the trees. This puzzling observation 
can be understood by comparing the relaxation time of the perturbations of pa{X, Y) along Y, Td ~ A^/D, with the 
time needed by the front of copper to move on the same length, r/ '-^ X/Xf = \\fT jvsfD. Since t/ increases with 
time T, in the asymptotic regime it will be much greater than r^. Therefore, whereas ps is highly correlated along the 
Y direction due to the structure of the solid deposit, there are eventually no fluctuations in pA along this direction. 

4. Departure from pure diffusion in the reaction zone 

The fact that the transference number i+ deduced from t+ = 1 — g and the inferred value oi q = F{v) from Eq. (p9) 
decreases significantly from 0.33 to 0.21 when the concentration of CuCl2 is increased from 0.5 mol.l"^ to 1 mol.r^^ 
is unlikely to be caused solely by a pure salt-concentration effect. It is also possible that convection produced by 
the sedimentation of CuCl crystallites toward the bottom of the cell could artificially increase the effective diffusion 
coefficient close to the reaction front by convective mixing. This would cause an increase of v{q) (the prefactor for the 
speed of the reaction front) which could at least partly explain the difference in the inferred q values, and therefore 
also in the effective i+ values. 

CONCLUSION 

We have shown that after long times the corrosion of highly porous copper clusters can be understood as a one- 
dimensional, homogeneous, mean-field A-l-B — > C reaction-diffusion process with one diffusing and one static reactant. 
This is the first experimental analysis of such a situation where only one reactant is free to diffuse through the other 
one. Whereas one would expect highly complex dynamics and a possible breakdown of the mean-field approximation 
when the reaction is confined to a porous (fractal) interface, we show that in this particular corrosion system, the 
dynamics are equivalent to those expected for an homogeneous system. The strength of our demonstration is built 
on precise measurements of the concentration field of the diffusing species by interferometry which are compared 
quantitatively with analytical predictions of the one-dimensional mean-field model. 
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FIG. 1. Schematic diagram of the thin-gap electrodeposition cell containing a ramified, metallic copper deposit. Note that 
the size of the deposit has been enlarged for clarity. 
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FIG. 2. (a) Photograph of a copper deposit grown from a 0.5 moLP CuCb solution at j = 40 mA.cm^ for approximately 
15 minutes, (b) Photograph of the same deposit half an hour after the current had been switched off. (The white zone is CuCl.) 
(c) The montage shows a sequence of photographs of a small region of the deposit including the reaction front taken every 30 
minutes after the interruption of the current. 
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FIG. 3. Log-log plot of the position of the reaction front Xf 
in CuCl2 0.5 mol.l"^ for deposits grown at j — 40 mA.cm"'^. 
one- dimensional, mean-field theory, given by Eq. ([ll|), with D = 



as a function of time T for two different experimental runs 
The solid lines of slope i represent the predictions of the 
10 cm .s , in the cases q = 0.6 and q = 0.73. 
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FIG. 4. Interferometric characterization of the concentration field around a copper deposit during its dissolution (a) just 
before the interruption of the current, (b) 15 minutes and (c) one hour later. (Ap « Pa/10 between adjacent isoconcentration 
lines.) The deposit grown in 0.5 mol.P^ CuCb solution at j = 40 rnA.cm"^ for 20 minutes. The concentration of CuCb is 
negligibly small inside and ahead (to the left) of the reaction front and approaches the bulk value of 0.5 mol.P^ far behind (to 
the right of) the front. 
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FIG. 5. One-dimensional concentration profiles extracted from the two-dimensional data. The deposit has been grown from 
a 1 moLP^ CuCl2 solution at j = 68 mA.cm"^ during 15 minutes. The concentration profiles are shown every 15 minutes after 
the current had been switched off. The different symbols are added on each profile to differentiate the recording times. These 
symbols will be used on the next representations of the concentration profiles in Figs 9 and 10. 
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FIG. 6. Log-log plot of the temporal evolution of the derivative of a{X, T) at X = as a function of T. Same parameters as 
in Fig. g. The plain line corresponds to the prediction of Eq. ( |4£| ) with D = 10~^ cm'^.s"^ and q = 0.79. 
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FIG. 7. The asymptotic similarity function a{x,i) ~ ^(C) where C = ^^ + v shown for Ao = and v = Q, 0.5, 1, 1.5, 2, 2.5 
from left to right. 
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FIG. 8. The exact asymptotic dependence of v, the square root of the dimensionless diffusion constant of the reaction front, 
on the asymmetry parameter q, predicted by Eq. (pSJ). 
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FIG. 9. Collapse of the experimental concentration data in the diffusion layer plotted versus the similarity variable 
{X — Xf)/2\/DT compared with the theoretically predicted asymptotic experimental similarity function A{Q (the solid line). 
The profiles are the same as those plotted in Fig. H, but only one point out of 20 is shown on this plot for clarity. 
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FIG. 10. Collapse of the experimental concentration data in the reaction front plotted versus the similarity variable 
{X — Xf)/2. The solid line shows the linearized extension of the similarity function A(C) from the diffusion layer (see Fig. 0) 
extended into the reaction front. These profiles are the same as those plotted in Fig. n, but only one point out of 4 is shown 
on this plot for clarity. The negative concentration values are artifacts of the interferometric technique and have no physical 
meaning. 
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FIG. 11. Relaxation of the two dimensional initial concentration field at the beginning of the dissolution, (a) Isoconcentration 
line a = 0.1, for T = 27 s. The deposit has been grown from a 1 mol.l"^ CuCl2 solution, at j = 68mA. cm~^ during 15 minutes, 
(b) Log-linear plot of the evolution of the amplitude G of the modulation of A concentration, as shown in (a), versus the 
reduced time T/rf = 4u'^DT/G{0f. 
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